Absolute Drell-Yan Dimuon Cross Sections in 800 GeV/c pp and pd Collisions 
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The Fermilab E866/NuSea Collaboration has measured the Drell-Yan dimuon cross sections in 
800 GeV/c pp and pd collisions. This represents the first measurement of the Drell-Yan cross section 
in pp collisions over a broad kinematic region and the most extensive study to date of the Drell-Yan 
cross section in pd collisions. The results indicate that recent global parton distribution fits provide a 
good description of the light antiquark sea in the nucleon over the Bjorken-a; range 0.03 < x < 0.f5, 
but overestimate the valence quark distributions as x — » f . 

PACS numbers: 13.85.Qk, 14.20.Dh, 12.38.Qk, 24.85. +p 



Drell-Yan dimuon production provides valuable infor- 
mation about the partonic structure of hadrons that is 
complementary to deep-inelastic scattering (DIS) stud- 
ies because it distinguishes between quarks and anti- 
quarks. The Drell-Yan cross section for pA^ fi + /i~X 
may be written in terms of the parton distribution func- 
tions (PDFs) of the colliding hadrons as 



dM dxF 9 x\ + X2 

X '^2 e q[n( X l)Q2(x2)+qi(xi)q2(x 2 )} 

1 

where the subscript 1 (2) denotes the beam (target) 
hadron. Measurement of the invariant mass, M 2 — 
X1X2S, and Feynman-x, xp — 2pi,/^/s — xi — £2, of 
the muon pair (where pl is the longitudinal momentum 
of the muon pair, and yfs the center-of-mass energy of 
the hadrons) determines the momentum fraction £1(2) of 
the beam (target) parton. Equation is exact to lead- 
ing order, but its general features are also preserved by 
next-to-leading-order (NLO) calculations. 

The Fermilab E866/NuSea Collaboration has mea- 
sured the Drell-Yan cross sections in 800 GeV/c pp and 
pd collisions. This represents the first study of the Drell- 
Yan cross section in pp collisions over a broad kinematic 
range and the most comprehensive measurement in pd 
collisions. Previous publications 0, Q have described 



the E866 measurement of the d/u ratio in the proton as 
a function of x, based on the X2 dependence of the Drell- 
Yan cross section ratio <r pd /2a pp , and those results have 
been included in recent global PDF fits 0, 0, H, H, . 

We present measurements of the pp and pd Drell- 
Yan dimuon cross sections for pair mass in the ranges 
4.2 < M < 8.7 GeV or 10.85 < M < 16.85 GeV and 
—0.05 < xp < 0.8. At the large values of Xf appropri- 
ate for much of the E866 data, the first term in the sum 
in Eq. ^ dominates and the cross sections are primarily 
sensitive to the valence distributions in the proton beam 
and the antiquarks at small x in the proton and deuteron 
targets. Thus, these measurements probe the partonic 
structure of the nucleon in two important kinematic re- 
gions. Recent global PDF fits use the precise HERA DIS 
cross sections to fix the magnitude of the light quark sea 
at small x and the E605 pCu cross sections |8| to fix 
the antiquark distributions for 0.14 < x < 0.3. In con- 
trast, Ref. |9j noted that the magnitude of the antiquark 
distributions near x ~ 0.04 is relatively unconstrained 
by current data, and this limits the precision of calcula- 
tions of the W production cross section at the Tevatron. 
There has also been considerable recent interest in under- 
standing the valence quark distributions in the nucleon 
as x — ► 1, which play a key role in searches for physics 
beyond the Standard Model, e.g., production of addi- 
tional vector bosons as predicted by left-right symmetric 
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models. To date, this region is constrained only by DIS 
cross sections, some of which suffer from limited statis- 
tics, while many others involve substantial ambiguities 
associated with corrections for the nuclear targets used 
0>0 (See Ref. for a recent review.) The present 
results provide new constraints on the magnitudes of the 
antiquark sea for 0.03 < x < 0.15 and of the valence 
quarks for x — ► 1. 

E866 used a 3-dipole magnetic spectrometer employed 
previously in E605 @, E772 and E789, modified by the 
addition of new detectors at the first tracking station. 
An 800 GeV/c proton beam bombarded identical target 
flasks containing liquid hydrogen, liquid deuterium and 
vacuum that were alternated every few minutes. After 
passing through the target, the beam was intercepted by 
a copper beam dump that was followed by a thick hadron 
absorber, ensuring that only muons traversed the spec- 
trometer's detectors. Drell-Yan events were recorded us- 
ing three different spectrometer magnet settings, chosen 
to focus low-, intermediate- and high-mass muon pairs 
and provide acceptance from below the Jji\) to above 15 
GeV. A detailed description may be found in Ref. 0- 

The present analysis used 55,000 pp and 121,000 pd 
Drell-Yan events, approximately half the statistics of the 
E866 d/u study |2|. The d/u analysis was optimized 
to achieve minimum relative uncertainties in a p /2a pp , 
whereas the present analysis was optimized for absolute 
measurements of the cross sections. Therefore, more 
stringent fiducial cuts were adopted in the present analy- 
sis, eliminating events for which the absolute acceptance 
of the spectrometer could not be reliably determined. 
This minimized systematic uncertainties at the cost of 
statistical precision, especially for X2 ^ 0.15. In con- 
trast, 7% of the events accepted in this analysis came 
from data sets that were excluded previously 0,0 . Sev- 
eral improvements were made to the event reconstruc- 
tion, notably involving the treatments of energy loss and 
multiple scattering in the absorber, that led to a more 
precise reconstruction of the dimuon kinematics and a 
better match between the real and simulated events. 

The spectrometer acceptance was calculated sepa- 
rately for each data set as a function of mass, xf, and 
transverse momentum (pt) using a detailed Monte Carlo 
simulation of the spectrometer. The virtual photon az- 
imuthal production and decay angles were assumed to be 
distributed isotropically, and the polar decay angles were 
assumed to be distributed according to 1+cos 2 9d, consis- 
tent with theoretical expectations and previous Drell-Yan 
angular distribution studies 0,0,0. Drell-Yan Monte 
Carlo events were thrown with realistic kinematic distri- 
butions in mass, xf, and pr, then reweighted to provide a 
precise match to the observed experimental distributions. 
Triply-differential cross sections, d 3 a/dMdxFdpT, were 
calculated for each kinematic bin, then integrated overp-r 
to obtain the invariant cross sections, M 3 d 2 a/dMdxF- 
The transverse momentum acceptance of the spectrome- 



ter extended to px ~ 5 — 7 GeV/c, but it was nonetheless 
necessary to extrapolate the integration to large px where 
the Monte Carlo showed zero acceptance. The extrapo- 
lation contributed well under 1% in most cases, and was 
never more than 5%. 

A detailed study was performed of the point-to-point 
systematic uncertainties in the measured cross sections. 
The two dominant contributions were the statistical un- 
certainty in the Monte Carlo event samples that were 
used for the acceptance calculations and the absolute 
field strength of the spectrometer magnets. Smaller con- 
tributions came from the uncertainties in the hodoscope, 
drift chamber, and trigger efficiencies, the composition 
and density of the targets, and the extrapolations to large 
Pf- Finally, some large- x p events contained a muon that 
passed very close to the edge of the beam dump, increas- 
ing the uncertainty in the energy loss and multiple scat- 
tering corrections. An additional ± 5% systematic un- 
certainty was added to the final results for the affected 
kinematic bins. The total point-to-point systematic un- 
certainties within any (M,xf) bin were strongly corre- 
lated between the pp and pd cross sections. In addi- 
tion to the point-to-point systematic uncertainties, there 
was a ± 6.5% overall normalization uncertainty, associ- 
ated with the calibration of the beam intensity. See Ref. 
0for additional details. 

Figure^shows the Drell-Yan pp and pd invariant cross 
sections per nucleon for selected Xf bins. The results 
agree with previous 800 GeV/c Drell-Yan cross section 
measurements in pCu collisions by E605 ||. They also 
agree with previous measurements in pd collisions by 
E772 for x F < 0.3. At larger x F and small M, 
the E772 cross sections are systematically larger than 
the present results. E772 quoted a larger point-to-point 
systematic uncertainty in this kinematic region, but the 
results of the two experiments differ by more than the 
combined systematic uncertainties would predict. The 
largest differences are in the region 0.4 < xf < 0.6 and 
4.2 < M < 7 GeV, where inconsistencies between the 
E772 results and expectations from global PDF fits were 
noted previously || . This region was studied during E866 
with very different acceptances by the low- and high-mass 
spectrometer settings, and the results are consistent. 

Figure^Jalso shows the results of next-to- leading-order 
calculations of the Drell-Yan cross sections based on the 
CTEQ6 and MRST2001 Q global PDF fits. The 
agreement with the global fits is very good over the entire 
kinematic region. This agreement may be quantified by 
computing a if' -factor, which we define to be the ratio 
of the experimental cross section to a NLO prediction. 
Table [I] shows i-T'-factors for several recent global PDF 
fits. With the exception of GRV98, all of the recent PDF 
fits predict the absolute magnitude of the Drell-Yan cross 
sections to within the ±6.5% normalization uncertainty. 

While the overall normalization is well reproduced, 
there are systematic deviations between the measure- 
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FIG. 1: FNAL E866 Drell-Yan cross sections per nucleon for selected if bins. The E866 pd (solid circles) and pp (solid 
diamonds) cross sections are shown in alternate decades, compared with previous pCu results from E605 |8J (open triangles) 
and pd results from E772 [l7j (open squares). NLO cross section calculations based on the CTEQ6 Q (dashed curves) and 
MRST2001 6] (solid curves) PDF fits are also shown. The error bars on the E605 and E772 data points are statistical only. 
Those on the E866 data points are the sum in quadrature of the statistical and point-to-point systematic uncertainties. An 
additional ± 6.5% normalization uncertainty is common to all E866 data points. 



TABLE I: /^'-factors obtained with various PDFs, where 
K' = a c * p /a NL °. The \ 2 values include the statistical and 
point-to-point systematic uncertainties in the experimental 
cross sections, but not the ± 6.5% global normalization un- 
certainty. The K' fits have 183 and 190 degrees of freedom 
(dof) for the pp and pd reactions, respectively. 



PDF 


K' 


X 2 /dof 




X 2 /dof 


CTEQ5 [31 


0.976 


1.42 


0.963 


2.51 


CTEQ6 [41 


1.016 


1.39 


1.001 


2.56 


MRST98 [5] 


0.973 


1.38 


0.960 


2.37 


MRST2001 [6j 


0.980 


1.45 


0.966 


2.44 


GRV98 [JJ 


0.811 


2.04 


0.808 


4.15 



merits and the predictions that are reflected in the large 
X 2 values. To elucidate these deviations, it is useful to ex- 
amine the experimental cross sections separately as func- 
tions of Xi and x 2 . Most of the events have x\ 3> x 2 , 
which implies that the x\ dependence is primarily sensi- 
tive to the valence quarks in the proton beam and the 
x 2 dependence to the antiquaries in the target. Fig- 
ure [5] shows the ratios of the experimental cross sec- 
tions to NLO calculations using the MRST2001 global 
PDF fit, separately as a function of x 2 and Xi, averaged 
over the other momentum fraction. The uncertainties in 
the NLO calculations from the PDF fit are also shown. 



After accounting for the normalization uncertainty, the 
MRST2001 partons provide a very good description of 
the X2 dependence of the pp cross sections over the full 
range, and a good description of the pd cross sections for 
x 2 < 0.15. The CTEQ6 fits describe the x 2 dependence 
equally well. This indicates that the current PDFs inter- 
polate the u and d distributions successfully between the 
HERA measurements at small x and the E605 measure- 
ments for x > 0.14, but these new results provide tighter 
constraints on the magnitude of the antiquark distribu- 
tions for 0.03 < x < 0.15 than have existed to date 0. 

Figure^ shows the ratios of the cross sections to NLO 
calculations using the MRST2001 PDF fits as a function 
of x\. The qualitative behavior of the CTEQ6 distri- 
butions is quite similar: both PDF fits overestimate the 
valence quarks by 15-20% at large Xi, and the dis- 
crepancy between the data and current PDFs appears 
to be larger for the pd cross sections than for the pp 
cross sections. The acceptance of these data includes an 
anti-correlation between x\ and x 2 0, so most of the 
events at large x\ also have x 2 < 0.1. The NLO calcu- 
lations describe the average x 2 dependence of the cross 
sections quite well for x 2 < 0.1, so the discrepancy in 
Fig. [2t> is isolated to the subset of the I0W-X2 events 
that also have large x\. Further comparisons between 
the measured cross sections and the NLO calculations, 
in which the x± and x 2 dependences are examined for 
limited ranges of the other variable [l3j. reinforce this 



4 



1.25 



0.75 



1.25 



D 

T3 



0.75 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- a) 


_ 


^^^^^^^^ 
* 

- +/- 6.5% normalization uncertainty 

-, , 1 , , , , 1 , , , , 1 , , , , 1 , , , , 1 




0.05 0.1 0.15 0.2 0.25 

X 2 


| i i i i | i i i i | i i i i | i i i i | i i i i | 

MRST2001 J 
' uncertainty 

pd fit 

— ppftt 


— * T'--i6— -4-.-- £ ^— _ 

* " » i * 

- • pd—>\l\L X """""*"--•-.. 

- A pp->\?)lX 

- +/- 6.5% normalization uncertainty 

T I T . T . T , T , T 1 , T T T ( T , T , T , 


, , 1 T 



0.3 



0.4 



0.5 



0.6 



0.7 



0.8 



FIG. 2: Ratios of the measured Drell-Yan pp (open trian- 
gles) and pd (solid circles) cross sections to NLO calculations 
based on the MRST2001 |t| PDF fit plotted vs. a) x 2 and 
b) xi, averaged over the other variable. The inner error bars 
represent the statistical uncertainty in the ratio, the outer 
error bar the sum in quadrature of the statistical and point- 
to-point systematic uncertainties. Solid lines represent the 
experimental uncertainty ranges 9] ona) u + J and b) 4u ± d 
in the MRST2001 PDF fit. Dashed and dot-dashed lines are 
the pd and pp fits described in the text. 



TABLE II: Fits of data/theory to the phenomenological form 
K'(xi) = a (1 — xi) 13 . Uncertainties on the fits include only 
the statistical uncertainties in the data. 





MRST2001 


CTEQ6 




pp pd 


pp pd 


a 


1.06 ± 0.02 1.09 ± 0.01 


1.11 ±0.02 1.15 ±0.01 


P 


0.14 ±0.03 0.21 ±0.02 


0.15 ± 0.03 0.23 ±0.02 



conclusion. Table [H] shows the results of fits to the X\ 
dependence of the i^'-factors with a phenomenological 
form, K'{xi) = a (1 - Xl f, for both PDF sets. 

Unlike large- a; DIS off deuterium targets, these data 
are in a kinematic region where nuclear-dependence cor- 
rections are known to be small [l8j . The pp and pd cross 
sections at large x\ constrain slightly different linear com- 
binations of u{x) and d(x), with greater sensitivity to u 
quarks. The results imply that the u quark distribu- 
tions in CTEQ6 and MRST2001 are overestimated as 
x — ► 1. Recent fits to HI charged- and neutral-current 
data have also indicated that the global PDF fits appear 
to overestimate the u quark distribution near x ~ 0.65, 
albeit with limited statistics 0]. The Drell-Yan cross 
sections may also point to problems with the d/u ratio 



as x — > 1 [ToL Hlj. However, determination of the d/u 
ratio at large x from the present data will require a fit to 
the full two-dimensional invariant cross sections together 
with the rest of the current world data. 

Figure[2]shows that the present results fall well outside 
the experimental uncertainty 9] for the valence quarks 
in the MRST2001 PDF fit. A similar comparison shows 
that the present results track the lower limit on the ex- 
perimental uncertainty 3 in the CTEQ6 PDF fit. These 
discrepancies reflect additional uncertainties in the PDF 
fits [lEl, and they imply that future PDF fits will see a 
substantial correction to the u and d quark distributions 
at large x. The reduced valence quark distributions at 
large x implied by these results will also lead to predic- 
tions of smaller cross sections for some signatures for new 
physics that are being sought at the Tevatron. 
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